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Abstract 

In this article, we estimate hydrostatic stress developed in gallium ion implanted 
gallium nitride epitaxial layers using Raman measurements. We have calculated 
deformation potential constants for i?2(high) mode in these epi-layers. The pres- 
ence of a polar phonon-plasmon coupling in these systems has also been demon- 
strated. In as-implanted samples, with an increase in implantation fluence, we 
have observed disorder-activated Raman scattering. 
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1. Introduction 

Wurtzite structure Gallium Nitride (/i-GaN) epitaxial layers (epi-layers) are 
direct wide band gap (~ 3.4 eV) two dimensional semiconductor systems which 
have applications in optoelectronic and microelectronic devices operating in the 
blue and ultraviolet. In GaN-based devices, ion-implantation is an attractive 
technique for selective-area doping, precise control over dopant concentration 
etc. Above a threshold fluence, ion-implantation leads to formation of cubic 
(c-GaN, space group T^) phase, which has many advantages over /i-GaN (space 
group Cq v ) £Q. The ion irradiation knocks out atoms from the irradiated ma- 
terial and creates point defects in the lattice. It does not create any massive 
damage and large defect complexes. It is to be noted that the important point 
defects, caused by ion-irradiation, are vacancies and interstitials. These defects 
induce localized electron energy levels into the band gap of the host semicon- 
ductor and hence results in a change in the electrical properties of the material. 
In addition, these defect states interact with light, inducing an increase in ab- 
sorption or emission of photons in radiative recombination processes 0J 
It has been claimed that the atomic structure of the defect states are respon- 
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sible for the yellow emission in GaN Strain in the lattice (due to lattice 
mismatch, difference in thermal coefficient, incorporated impurities, and point 
defects) is an important issue in fabricating optoelectronic devices. 

Stress-strain relation in virgin /i-GaN layer is well-studied in the literature 
7 . In this article, we discuss Raman measurements on self- ion (Ga ++ ) im- 
planted GaN epi-layers, with an emphasis on induced stress in the system. From 
our experimental results, we have estimated the hydrostatic strain coefficients 
inside the layers. In addition, we show the effect of disorder-activated Raman 
scattering in as-implanted samples and demonstrate the presence of a polar 
phonon-plasmon coupling in post-annealed implanted epi-layers. 
2. Experimental Details 

GaN epi-layers (pristine sample), of thickness 6 /im, were grown on (0001) 
AI2O3 substrate in a horizontal Metalorganic Chemical Vapor Deposition reac- 
tor at atmospheric pressure. Trimethylgallium (TMGa) was used as precursor 
with NH3 as reactant gas. The flow rate of NH3 was 3000 seem. The low- 
temperature-deposited GaN buffer layer was first grown at 550C for about 30 
min. Temperature was raised to 1050C and growth is continued for 8-10 hrs. 
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Ga ++ ions were implanted in these layers at 3.0 MeV with different fluences, 
1 x 10 15 , 2 x 10 15 and 1 x 10 16 per cm 2 (Sample A - Sample C). The damage of 
the crystal lattice, induced by implantation, was removed by thermal annealing 
in flowing ultra-high pure N2 at 650 °C for 15 minutes and subsequently at 
1000 °C for 2 minutes The undoped and implanted epi-layers were found to 
be of n-type from Hall measurements with carrier concentration ~ 4-8 x 10 17 
cm~ 3 . Raman measurements were carried out at room temperature in a back- 
scattering geometry using a 488 nm air cooled Ar + laser as an excitation source. 
To avoid phase change in the samples by laser heating, the power density on 
the samples was tuned to 3 x 10 4 Watt/m 2 . The Raman spectra were obtained 
using TRIAX550 single monochromator with an open electrode charge coupled 
device as a detector. With 150 fim slit-width of the spectrometer, the accuracy 
of our Raman measurement is ±1.5 cm -1 . 
3. Results and Discussion 

3.1. Pristine and self-ion implanted GaN epi-layer 

Fig. 1 shows the Raman spectra of pristine and as-implanted GaN layers 
over a wavenumber range of 240 to 850 cm . 

4 



Pristine GaN epi-layer: 

Group theory predicts that for h-GaN there are eight phonon modes, 2Ai + 
2£>i + 2Ei + 2E2, at the T point. Except one A\ and one E\ mode (which 
are acoustic), the remaining six modes are optical. Optical A\ and E\ modes 
are both Raman and infrared active, while the two Ei modes are only Raman 
active and the two B\ modes are silent (Raman inactive). The higher frequency 
branch of the E2 phonon mode is denoted by E^ [H|- Peaks at 374, 416 and 
751 cm" 1 , shown by * marks in Fig. 1, are phonons from the sapphire substrate 
|1(J|. In the Raman spectrum of the pristine sample the line at 568 cm -1 is due 
to the E2 mode. The Raman mode at 345 cm -1 is not allowed by the C6„ 
space group in the first order Raman scattering at the zone center. This peak 
can be attributed to the 2nd order Raman scattering due to acoustic overtones 
in this region (300 - 420 cm^ 1 ) |11|. In the high frequency region for the pristine 
sample, we distinctly observe a feature at 743 cm -1 (shown by double arrows 
in Fig. 1) due to the polar mode. This appears as a shoulder of the prominent 
peak at 751 cm -1 from sapphire. The origin of this peak will be discussed latter. 
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Self-ion implanted GaN epi-layers: 

After ion-implantation the system loses its crystallinity and becomes disordered. 
This is reflected in the broadening of the spectra of ion-implanted GaN layers 
(Fig. 1). In addition, for as-implanted samples, we observe a broad band 
at around 290 cm -1 , which becomes prominent at higher fluences (inset of 
Fig. 1). As a result of ion- irradiation, there is a slight rearrangement in the 
lattice structure, hence, the wavevector conservation rule in Raman scattering 
gets relaxed. The Raman spectrum reflects the total phonon density of states 
(so called disorder-activated Raman scattering). The Raman-inactive modes 
become Raman-active. Though silent for bulk h-GaN, the phonon density of 
states of the highest acoustic phonon branch at the zone boundary, the B\ 
mode, has a strong feature in the region 290-300 cm -1 ^J^HEl- We assign 
the above observed feature at around 290 cm -1 to this B\ mode pp. Due to 
disorder-activated Raman scattering, the B^ mode also becomes allowed and 
appears as a small hump around 674 cm -1 , shown by o mark in Fig. 1 (See 
phonon density of states in Ref. The polar mode at around 740 cm -1 is 

also broadened. In these as-implanted samples, the absence of Raman lines from 
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sapphire (substrate) indicates the effective "trapping" of the laser light inside 
the layers, which can affect the Raman scattering from the layers by changing 
the selection rules. 

3.2. Post-annealed implanted samples 

Raman spectra of the post-annealed GaN layers are shown in Fig. 2. After 
annealing, the disorder in samples gets reduced. All Raman modes in the low 
frequency region, discussed for the pristine sample, reappear distinctly in their 
spectra. In Raman spectra of post-annealed implanted samples the line at 573 
cm -1 is due to E% mode, shifted by 5 cm^ 1 from its value (568 cm -1 ) in bulk 
|JJ. However, we do not observe any variation in spectral line-shape and position 
with implantation (Sample A - Sample C). Peaks from the sapphire substrate 
clearly reappear in these spectra. With the decrease in disorder in the annealed 
samples, the intensity of the peak at 290 cm -1 due to the disorder-activated 
Raman scattering also decreases (inset of Fig. 2). However, we would like to 
point out that for c-GaN LA(X) and LA(L) (LA: longitudinal acoustic) modes 
appear at 286 cm -1 and 296 cm -1 ^3]. The photoluminescence spectra of our 
samples (not shown here) do not carry a clear signature of the presence of cubic 
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phase in them. 

3.3 LO phonon -plasmon coupling 

In Fig. 3 we have shown Raman spectra in [z(y,y)z] and [z(x,y)z] configurations 
for Sample A. The low frequency acoustic mode at 345 cm^ 1 and high frequency 
polar mode at 743 cm -1 appear only in parallel polarization. Similar polariza- 
tion dependence of the polar modes has been observed for all post-annealed 
implanted samples. In h-GaN, polar Ai(LO) [LOdongitudinal optical] mode is 
expected to appear at 733 cm -1 in [z(y,y)z] configuration only. Assigning the 
observed feature at 743 cm -1 to this Ai(LO) mode is not justified, because the 
shift of a Raman line by 10 cm -1 cannot be explained either by strain or zone 
folding effect [T2|. Also, it is to be noted that this is the only peak for which 
the frequency shift is more than the shift for the other peaks in Fig. 2. 

Undoped or nominally doped GaN is invariably of n-type, usually with a 
high free electron concentration (10 17 -10 18 cm~ 3 ) at room temperature |17|. 
The coupling of collective electronic excitation (plasmon) and lattice vibration 
gives rise to the LO phonon-plasmon coupled mode. Such a mode in h-GaN 
shifts the LO phonon to the higher frequency HJ. Assigning the Raman line at 
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743 cm 1 to the plasmon-Ai(LO) phonon mode, we estimate the charge carrier 
density in Sample A using the empirical relation 

n = 1.1 x 10 17 (Ao;) - 764 , (1) 

where, Au> is the observed frequency shift of the LO phonon. Taking Alj = 10 
cm -1 , the n is estimated to be 6 x 10 17 cm~ 3 . which is also the value obtained 
from the Hall measurement. The above equation is valid for n < 1 x 10 19 cm -3 . 

Considering the contribution of the deformation potential within the lattice, 
a more precise evaluation of n can been carried out. The intensity of the Raman 
scattering is then given by [HI ED D2| 

I{u)) = SA{uj)lm[-l/e{uj)l (2) 

where, lu is the Raman shift and S is the proportionality constant. The dielectric 
function is 

e(uu) — [1 + [lo\ — Wy)/ (wj — uj 2 — iuoT) — lu 2 p / (oj 2 — itiry)] = 0. (3) 

Here, lot and wl are the transverse optical (TO) and LO phonon frequencies for 
the uncoupled phonon modes, eoo is the high-frequency dielectric constant of 
GaN. r(7) is the damping rate of the phonon (plasmon). The plasmon frequency 
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up and A(uj) are given by 

lu p = [A^/ie^m*)] 1 ' 2 (4) 

and 

A(lo) = 1 + 2Cc4[u^7(<4 - uj 2 ) - lu 2 T{lu 2 + 7 2 - w|,)]/Aw 
+C 2 (4/Aw)[^[ 7 K - uj 2 t ) + I>f, - 2lu 2 )} + lu 2 T(lu 2 + 7 2 )]/(w| - (5) 

Here, C is the Faust-Henry coefficient, can be estimated from the ratio of the 
intensities of the polar LO and TO Raman scattering modes measured at 90° 
geometry [2D] ■ The effective mass of the electron is m* = 0.2mo. All three 
parameters, T, 7 and n, are free parameters, which need to be chosen correctly, 
by fitting more than one Raman spectra of the layers with different carrier 
concentrations. Unfortunately, the Raman line shapes of all our post-annealed 
samples are nearly same, except, with higher implantation dose, a slight decrease 
in intensity of the coupled phonon-plasmon mode (modified A\(LO) mode) is 
observed due to damping by the hole plasmon. In such a case, it was not 
possible to choose all three parameters , T, 7 and n meaningfully, to estimate 
more accurate value of n by following this procedure, compared to what we have 
obtained from Eqn. 1. 
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3.4 Stress Analysis 

An epi-layer, of thickness 6 /jm, is not expected to exhibit a shift in the 
Raman line due to the effect of confinement of phonons |2J ■ Thus, the origin of 
the shift in Raman frequency from bulk samples can be either due to the presence 
of elastic strain or due to the change in interaction between the elastic medium 
and the macroscopic field. To eliminate the second factor we have investigated 
the non-polar E% mode. It is important to note that within our experimental 
accuracy, all annealed Ga ++ implanted GaN epi-layer exhibit nearly identical 
Raman spectra. Thus, the following stress analysis holds good for all post- 
annealed implanted samples. 

GaN epi-layers grown on sapphire contain residual strains produced by mis- 
match in the lattice constants and thermal expansion coefficients between the 
GaN film and substrate [22;. Sufficient residual strain in the film induces for- 
mation of dislocations and stacking faults. Residual strain also leads to wafer 
bowing. It has been shown in Ref. 23 that for the above defects in a GaN 
layer of thickness 1 (im (on sapphire substrate) the bi-axial strain energy varies 
from (GaN surface) to 0.4 GPa (GaN/sapphire interface), indicating that the 
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surface layer is free from this type of strain. On the other hand, the ion irradia- 
tion induces strain in the film due to point defects. These defects are expected 
to be more near the surface of the layer. As the thickness of our films is 6 fim, 
we attribute the shift in non-polar Raman mode only to the strain, generated 
by the point defects in the layer. These defects induce three dimensional stress, 
such as hydrostatic stress, if they are uniformly distributed. Based on the as- 
sumption of pure elastic theory (where Hooke's law is valid) the compressive 
stress can be calculated from the shift in E% from [21] 

LO — LUq = 4.17<7£f, (6) 

where, <jh is the hydrostatic stress in GPa, to and loq are the Raman line position 
of the stressed and unstressed samples. For a shift of 5 cm -1 of the E% mode 
between bulk GaN and ion-implanted layers we have estimated an ~ 1.20 GPa 
(see Table I). 

Knowing oh, the hydrostatic strain tensor components, u xx and u zz (inplanc 
and normal components of strain tensor) can be estimated from the following 
relations [23] 

U X X = Uyy = <7 H /Y (7) 
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and 

u zz = -Ru xx . (8) 

The Young's modulus, Y, of the material, in terms of elastic stiffness constants 
Cij is given by [21] 

Y _ (Cn + Ci2)C , 33 - 2Cf 3 

C33 — C13 

and i? is the Poisson's ratio, can be written as 



i?=^i±^-^. (10) 

L-33 — W3 



Using Eqn. 7 to Eqn. 10, the strain tensor components have been estimated 
and tabulated in Table I. Here, the elastic constants Cij are taken to be C\\ = 
396, C12 = 144, C13 = 100, and C 33 = 392 GPa [13 . 

The change in frequency for a given phonon mode under symmetry conserv- 
ing stress can be expressed in terms of strain tensor components and deformation 
potentials constants (a a and 6a) as [22 

oj - Uq = 2a x u xx + b x u zz . (11) 

Furthermore, the bulk Grueneisen parameter 7a is related to the characteris- 
tic phonon frequency at hydrostatic compression, phonon deformation potential 
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constants, and elastic constants of ft.-GaN by [7j 



7A = - 



2a x {C 33 - C 13 ) + 6a(Cii + C12 - 2C 13 ) 
co x (C 11 + C 12 +2C 33 -4C 13 ) 



(12) 



For GaN, C 33 + C\ 3 = C\\ + C\ 2 . Therefore, from the above relation we get, 



The Griieneisen parameter used is j e h = 1.54 7 . Knowing the values of strain 
tensor components, the deformation potential constants have been calculated 
to be -81 cm- 1 and -2486 cm" 1 . 

With irradiation, we do not observe any noticeable change in the peak po- 
sition and lincshapc of non-polar Raman modes, which implies that no extra 
stress has been introduced in layers with increase in implantation fluence. The 
values of above parameters u xx , u zz , a\ and b\ are nearly same for all implanted 
samples. 

3.5 Conclusion 

We have performed first-order Raman measurements on Ga ++ implanted GaN 
epi-layers. Raman lines become broadened after implantation. However, they 
reappear in all post-annealed implanted samples. We do not observe any change 



7a = 



2aa + b x 
3uj\ 



(13) 
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in peak positions in Raman spectra of the films with an increase in implantation 
dose. Thus, the estimated hydrostatic stress and strain tensor components, 
obtained by us from the knowledge of phonon modes, are expected to hold 
good for all post-annealed implanted samples. It would have been interesting 
to check the stress analysis by high resolution x-ray diffraction measurements. 
We have also the disorder-activated Raman scattering in as-implanted samples. 
Interesting behavior of polar modes is reported for the post-annealed implanted 
samples. 
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Figure Captions 

Figure 1. Raman spectra of pristine (P) and as- implanted GaN epi-layers (Sam- 
ple A - Sample C). The region between 240-340 cm -1 is shown in the inset. 

Figure 2. Raman spectra of pristine and post- annealed Ga ++ implanted GaN 
epi-layers (Sample A- Sample C). The region between 240-340 cm" 1 is shown 
in the inset. 

Figure 3. Polarization dependence of Raman modes in Sample A. Other im- 
planted samples exhibit similar behavior. 
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Table Caption 



Tabic I. Hydrostatic stress, phonon deformation potential constants and strain 
tensor components of postannealed Ga ++ implanted GaN epi-layers. 
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Phonon 






ax 


bx 




u zz 


Symmetry 


(cm -1 ) 


GPa 


cm -1 


cm -1 


x (10~ 3 ) 


x (10~ 3 ) 


£ 2 (highj 


5 ±1.5 


1.20 ±0.36 


-81 ±1 


-2486±1 


1.83 ±0.55 


-2.13 ±0.64 



Table 1: Hydrostatic stress, phonon deformation potential constants and strain 
tensor components of Ga ++ implanted GaN epi-layers. 
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Figure 1. S. Mai etal. 




Figure 2. S. Mai et al. 
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Figure 3. Mai et al. 



